Optimal quantum searching to find a common element of two sets 
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Given two sets A and B and two oracles Oa and Ob that can identify the elements of these sets 
respectively, the goal is to find an element common to both sets using minimum number of oracle 
queries. Each application of either Oa or Ob is taken as a single oracle query. This is basically a 
search problem and a straightforward application of Grover's algorithm can solve this problem but 
its performance is slow compared to the optimal one by a constant factor of 1.57. Here we present 
a variant of Grover's algorithm which achieves the optimal performance in not too restrictive cases. 
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I. INTRODUCTION 

In a general setting, Grover's algorithm [l| drives a 
quantum system from a known initial state \s) to an un- 
known final state which are also called source and 
target states respectively. This is done by iterating the 
operator G = Iglt on \s), where = 1 — 2\i(j) {ip] is the 
selective phase-inversion of the {ip) state. It is known that 
7r/4|(s|t)| iterations of G on \s) will bring the quantum 
system very close to the desired final state \t). Though 
the operator Ig is implemented using our knowledge of 
the source state |s), same is not true for It which has to 
be implemented using an oracle O that can identify and 
selectively apply a phase-inversion of the target state \t). 
Grover's algorithm is proved to be optimal 

In the case of search problem, we have a set of N items 
labeled by i € {0,1, . . . ,N — 1}. Each item corresponds 
to a basis state |i) of a quantum system having an A^- 
dimensional Hilbert space. We denote the set of solutions 
by t and let M = \t\he the total number of solutions. The 
source state is chosen to be the uniform superposition of 
all basis states \s) = (I/n/A") K)- The target state 
\t) is a uniform superposition of all target states \t) = 
(l/yM)E.etN)- So, \{s\t)\ = ^M/iVand ^^/nJM 
iterations of G are sufficient to get the state \t) which 
always gives a target state i d t upon measurement. To 
implement It, we define a binary function f{i) which is 
1 for i € t else 0. We attach an ancilla qubit, initially in 
state I a), to our search-space and then we use the oracle 
O to compute f{i) to get the transformation 
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Choosing \a) — |— ) = (|0) — \l))/^/2, above transforma- 
tion is equivalent to the It transformation on the search 

space, i.e. \i) ^ 

We consider the problem where we are given two sets 
A and B and two oracles Oa and Ob which can apply 
selective phase-inversions on the elements of sets A and 
B respectively. For this problem, it is natural to con- 
sider each application of either Oa or as a single 
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oracle query. Our goal is to find an element i e {AO B) 
which is common to both sets using minimum number of 
oracle queries. A quite related problem is the schedul- 
ing problem, where two parties Alice and Bob has a list 
of calendar dates on which they are free. They want to 
schedule their meeting by finding a date on which both of 
them are free. They want to do so by using least amount 
of communication. This problem was first considerd by 
Buhrman, Cleve and Wigderson (BCW), where they had 
shown that qantum mechanics can give a quadratic im- 
provement in the communication complexity of the prob- 
lem Q. Later, Grover came up with an algorithm with 
a better communication complexity [4]. 

Our problem is a simple variant of search problem 
and we can use Grover's algorithm to solve this prob- 
lem by implementing the It operator using the given or- 
acles Oa and Ob- As we show in next section, Grover's 
algorithm takes ^^JN/M oracle queries to solve this 
problem where M = |yl fl i?| is the total number of el- 
ements common to both sets. The optimality proof of 
Grover's algorithm can be used to pro ve tha t any al- 
gorithm will take at least gopt = 7r/4(-\/A^/M) queries. 
Because if an algorithm takes less than gopt queries then 
by choosing A — B, we can get a faster quantum search 
algorithm which is impossible. Since Grover's algorithm 
takes ^y^N/M = fgopt ~ l-57gopt queries for our prob- 
lem, its performance is suboptimal by a constant factor. 
In Sec. Ill, we present a variant of Grover's algorithm 
which achieves the optimal performance in not too re- 
strictive cases and then we conclude in Sec. IV. 



II. USING GROVER'S ALGORITHM 

Here we discuss straightforward applications of 
Grover's algorithm to our problem. We present two 
such algorithms based on different approaches. These 
algorithms are motivated by the quantum algorithms for 
scheduHng problem [3, 4]. 

Algorithm- 1 In our problem, < G [Af] B) and we need 
to implement It using given oracles Oa and Ob- Let us 
define two binary functions: fA{i) which is 1 for i £ A 
else and similarly, /_B(i) which is 1 for i G -B else 0. 
These functions can be computed using our oracles and 
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we have 



\^)\a®fA{^)) , \^)\a) ^\^)\a®fB{^)). (2) 



One possible way to implement It is to attach an ancilla 
qubit with the initial state \a) — |0). Then we use oracle 
Oa and apply Hadamard gate Hanc on the ancilla qubit 
to get 



z)|0) 



|0) + (-l)/-«|l) 



V2 



(3) 



The final state of ancilla qubit is | — 

erwise 1+) = (|0) + |1))/V2. Next, 

Ob- Since | + ®fB{i)) 
(„1)/b(«)U\ 



if fA{l) - 1 0th- 

we use the oracle 
= and I = 

we find that the basis state \i) gets a phase 



factor of -1 iff /^(i) = = 1, i.e. if i G {An B) 

is a target state. Then the ancilla qubit can be decou- 
pled from our system by doing the reverse of Eq. ([3]) and 
effectively we implement It on the search-space. This 
implementation uses 3 oracle queries, 2 of Oa and 1 of 
Ob, and since Grover's algorithm uses Qopt — j\/N/M 
applications of It , we find that a total of S^opt queries are 
needed to solve our problem using Algorithm- 1. 

Algorithm-2 A second approach is to first drive our 
quantum system from the uniform superposition of all 
basis states \s) — (I/Vn) \i) to the uniform superpo- 
sition of only those states corresponding to the elements 
of set A, i.e. \sa} — /! J2ieA N)- This is done by 
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treating i £ A states as effective target states in Grover's 
algorithm. The corresponding It transformation is then 
I A, a selective phase inversion oi i G A states, and it is 
implemented using the oracle Oa that computes in 
the same way as It was implemented in Eq. ([1]). We get 
the state \sa) using qa = n/^^/N/\A\) iterations of the 
operator Is I a on |s). So we have 
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Now let us consider \sa) state as an effective source state 
and iterate the operator Isa^b on \sa), where Ib is the 
selective phase-inversion of i G _B states, implemented 
using oracle Ob, and the operator Ig^ is given by 



IsA = I-2\sa){sa\ 

= I-2{IJaY^\s){s\{IaIsY^ 
= {IsIaY^Is{IaIsY^- 



(5) 



Since the state \sa) has close to zero amplitude in 
i ^ A states and Ib inverts the phases of only i G B 
states, we find that \{sA\t)\ = |AnB|/|A|, where \t) = 
{l/yy\A n B\)J2i£(AnB) K) target state. Iterat- 

ing IsA^B will amplify the amplitudes of only desired 
states and the number of iterations required to get |t) 
is 7r/4(-y/|A|/|ylni?|). Above equation implies that each 
iteration of Isa^b requires 2qA + 1 oracle queries, 2qA 



for /g^ and 1 for Ib- Hence the total number of oracle 
queries required by Algorithm-2 is 




Though this algorithm is better than Algorithm- 1, this 
is slow compared to the optimal performance by a factor 
of 7r/2 « 1.57. So far, we are trying to use Grover's 
algorithm to solve our problem. In the next section, we 
present a variant of Grover's algorithm which naturally 
relates to our problem and solves it optimally. 



III. VARIANT OF GROVER'S ALGORITHM 
A. The search operator 

Here we present a variant of Grover's algorithm for our 
problem. Our algorithm starts with the uniform super- 
position of all basis states \s) — {\/\/N) \i) and then 
iterates the following operator on it. 



V = UbUa- 



(7) 



To analyse the iteration of V on \s), we need to find 
its eigenspectra. We divide the set N of all basis states 
? S {0, 1, . . . , TV — 1} in four mutually non-intersecting 
sets: (1) Set A' whose elements are in set A but not in 
set B, i.e. A' = A - {An B), (2) Set B' whose elements 
are in set B but not in set A, i.e. B' = B — {An B), 
(3) Set T whose desired elements are common to both A 
and B, i.e. T = A n B, and (4) Set P whose elements 
are neither in set A, nor in set B, i.e. P = N — {AU B). 
Correspondingly, we define four quantum states \p), 
I a), and which are the uniform superposition of basis 
states of the corresponding sets, i.e. 



ieA' ieB' 

\t) = \T\-'/'j2\'), b) = |prV2^|,). 



ieT 



ieP 



These states are mutually orthogonal to each other and 
the initial state |s) can be written in terms of them as 

\s) ^a\a)+l3\b)+-/\t)+5\p), + +-f^ + 5^ ^ 1. (9) 

where 



m 

N 



(10) 



It is easy to see that our search operator V preserves 
the four-dimensional subspace of A^-dimensional Hilbert 
space. This subspace is orthogonally spanned by the 
states I a), |&), \t), and \p)- The operator I a inverts the 
phases of \a) and \t) and the operator Ib inverts the 
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phases of \b) and \t). The 4x4 matrix representation of 
these operators in this subspace are given by 



I A = 



-10 
10 
0-10 
1 



In = 



10 
0-100 
0-10 
1 



(11) 

To find the matrix representation of Is, we use the Eq. 
([9]) and the expression Is = I — 2\s){s\. We find it to be 



1 



- 2^2 -2al3 

-213a 1 - 2/32 

-27a -27/3 

-25a -2613 



-2q:7 
-2/37 
-272 
-2(57 



-2a5 
-2(35 
-27(5 
-2,52 



(12) 



Putting these expression in Eq. ([7]), we find that the 
matrix representation of our search operator V in the 
preserved four-dimensional subspace is given by 



V 



iPt - 1 4:a(3l2 
-AaPh 813'^h - 

4/3612 



1 



-4ck7Z2 
8a6r 



-4a7^2 

-8/37^1 
- 872^1 



where r = /32 + 72 and Ik — 1 
eigenspectra of V. 



—8a6r 
4/3612 
n 47(5^2 
-47^/2 1 - S6^r 

(13) 

~ kr. Next, wc find the 



B. Eigenvalues 

To find the eigenvalues of V, we note that this is a real 
operator, so its eigenvalues come in complex conjugate 
pairs. Further, we make an important observation that 



(14) 



The latter identity is because is a real unitary oper- 
ator, so V'^ — . Let be an eigenstate of V with 
eigenvalue e**, i.e. V\^p) 

{V-^V'')\t^) = (2cos0)|V') (15) 
(y-yt)l^) ^ (-2isin6')|V'). (16) 



e'^\^). Then 

(2cos0)|V') 
(-2isin6')|V'). 

Using Eq. ([TI)) . we get 

^[1 + {-ly+^v^J'lpJ ^2coBe^j 

Y^[l - (-l)''+^']V^jV'j = 2«sin6'V'j . 



(17) 



First of the above two equations can be written as eigen- 
value equations of two 2x2 sub-matrices, Vodd and Vi^von, 
of V as 



Vodd 

and 









Viz' 






— COS 6 








_ V^13 


^33 




_ ■(/'3 _ 


_ ^3 _ 






'V22 


V2a' 






= cos 9 








V2i 


Vii_ 




. "^4 _ 





(18) 



(19) 



It is easy to verify that the two eigenvalues, cos0±, of 
sub-matrices Vodd and V^von are same and they determine 
the four eigenvalues e^**± of V . We have 



cos a± 

X 

c 

?7 



(V^ii + V^33)/2 
("^"ll - ^33)74 
^3 • 



(20) 



Using expressions for {V\\, V13, V^33}, we get 

X = 4(a2 + 72)(/32 +^2)_4^2 ^ 

C = [l- 4(/32 + 72)(a2-72)_472 
T] = [4a7(2^2 ^ 272 - 1)]2 

Next, we find the eigenvectors of V 



(21) 



C. Eigenvectors 

The eigenvectors of V are found by finding the eigen- 
vectors of Vodd and T4vcn- Eqs. and gives the 
following relationships for the eigenvectors 



^"3 
^"4 



= ffiV'l , 9± = (COS0± - Fii)/Vi3 , 

= h^'^2 , h± = (C0S(?± - V22)IV2^ (22) 



Since the off-diagonal elements of V^dd and T4vcn are 
equal, we have g+g- — h^h- = — 1 so that we can write 
g± = ±5=*=^ and h± = ±/i^^ Also, Eq. ([T7]) implies that 
if we choose V'l to be real, then tp^ is also real while '02 
and -04 are imaginary. We also note that V\ip) = e^^jV') 
implies that V\ip*) — e~*^|V'*) since V is real. 

With all these observations, we can write down the di- 
agonalizer matrix D whose columns are the eigenvectors 
of F i.e.. 



D 



X X gx gx 

ly —ly ihy —ihy 

gx gx —X —X 

ihy —ihy —ly ly 



(23) 



and we have 



DVD'f 










(24) 



Here x and y are the arbitrary constants in the choice of 
eigenvectors for V^dd and Vcvcn respectively. To fix these 
constants, we use the normalization condition, i.e. choose 
the column vectors of D to be of unit length. Note that 
the first two column vectors of D, \Dii) and \Di2), are 
mutually orthogonal and the complex conjugates of each 
other. If we add them, then the length of sum \Dii) + 
\Di2) = 2{x\a) — x/g\t)) is 2. Similarly, if we subtract 
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them, then the length of the difference \Dii) — \Di2) 
2i{y\a) — y/h\t)) is 2. These give us 



1 



1 



; y 



where 



g = g+ = (cos 0+ - 1^11 )/Pi3 

h = h+ = (cOs6l+-y22)/V24 



(26) 



These equations determine the eigenvectors of V. Next, 
we analyse the performance of our algorithm. 



D. Performance 

To discuss the performance, we need to analyse the 
iteration of V on the initial state |s) given by Eq. 
This can be done by working in the eigenbasis of V. To 
write \s) in this basis, we multiply \s) by to get 



a\tp 



++/ 



K* lip- 



(27) 



where \4'±±) denote the normalised eigenvectors of V 
with the eigenvalues e**^± (for example ) is the nor- 
malised eigenvector with eigenvalue e"*"*^" ) and 



x{a + gj) — iy{/3 + hS) 
x{ad — 7) — ~ ^) 



(28) 



Thus we have 



= ae"i'+\il:++)+a*e-"i''+\i).+) (29) 
+ Ke*«^-|V'+-) + K*e-'«''-|V'— ). (30) 

Now, Eq. ^ gives 

\t) = gx{\ij++) + |^_+)) - x{\ij+^) + |^__)) , (31) 

and hence 

(i|T/9|s) = 2gx\(T\ cos(Zf7 + qe+) - 2x\k\ cos(Zk + qO-). 

(32) 

The probability of getting a target state after q iterations 
of V on the initial state is given by the square of above 
quantity. So far, we have presented the analysis of general 
case. Next, we consider the case when our algorithm 
achieves the optimal performance. 



After little calculation, we get 



(25) VCT^ = (l-4a2/32)+7' 



4(/32 - - 1) + 



Using last two equations in Eq. (f20l) . we get 



8a^(2/3^ - If 
1 - 4^2/32 
(35) 



cos^?+ = 1— 87 



1 - 4a2^2 



471 



l-a2_^2 



l-4a2^2 ' 
(36) 

and COS0- = — 1 + ^a^0^ + 0(7^). From now onwards, 
we ignore 0(a, /3, 7) terms in our calculations. Using Eq. 
(|26t. we get 



.9 = -(2a7)-i > 1; = 1/2/3. (37) 
Then Eq. ^ gives 



V2a-i < 1; y= ^2/3 . 



(38) 



and Eq. gives 



CT = -V%/2+0(/3^) =^ Zcr « -it/2: k = i/3/V2+0{l3^) . 

(39) 

Since 5 3> 1 and \cr\ > |k|, we find that the second term 
in the R.H.S. of Eq. (|32|) can be neglected and using 
above equations, we get 

=cos(|+4g7) +0(/3,a,7). (40) 

So we get a near to 1 amplitude in the target state after 
7r/87 iterations of our search operator V. Since each 
iteration of V takes 2 oracle queries, 1 for Oa and 1 for 
Ob, the total number of queries used by our algorithm 
is Qopt = which is also the optimal one. Thus, our 

algorithm is optimal as long as our problem parameters 
satisfy the condition psp which is, by no means, too 
restrictive. 



E. Optimal Performance 

We consider a special case when 

7<1, a,/3<l/10 (33) 

Note that 7 ^ 1 is not at all a restrictive case as if 7 is 
comparable to 1 , then we can easily find the desired item 
by random trials. As 7 ^ 1, we ignore 0(7'^) and higher 
terms in Eq. (I?T|) to get 

X = 4^2/32+ 472(^2 + /3^-l) 

ri = 16720^(2/32 - 1)2 (34) 



IV. CONCLUSION 

We have presented a variant of quantum search al- 
gorithm, which naturally relates to our goal of finding 
a common element of two sets. We have shown that 
our algorithm can be exactly analysed by considering a 
four-dimensional subspace and by using the properties of 
the matrix representing our search operator. Though a 
straightforward application of standard quantum search- 
ing will give a sub-optimal performance, our algorithm 
can give optimal performances in not too restrictive con- 
ditions. 
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